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Two-Band- Type Superconducting Instability in MgB2 
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Using the tight-binding method for the 7r-bands in MgB2 , the Hubbard on-site Coulomb inter- 
action on two inequivalent boron p^-orbitals is transformed into expressions in terms of vr-band 
operators. For scattering processes relevant to the problem in which a wave vector q is parallel 
to i, it is found to take a relatively simple form consisting of intra-band Coulomb scattering, 
interband pair scattering etc. with large constant coupling constants. This allows to get a simple 
expression for the amplitude of interband pair scattering between two 7r-bands, which diverges 
if the interband polarization function in it becomes large enough. The latter was approximately 
evaluated and found to be largely enhanced in the band structure in MgB2 . These results lead to 
a divergent interband pair scattering, meaning two-band-type superconducting instability with 
enhanced Tc- Adding a subsidiary BCS attractive interaction in each band into consideration, a 
semi-quantitative gap equation is given, and Tc and isotope exponent a are derived. The present 
instability is asserted to be the origin of high Tc in MgB2 . 

KEYWORDS: MgB2, interband polarization function, superconducting instability, two-band mechanism, interband 
pair transfer, Coulomb interaction, nesting, isotope effect 
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The recent finding oiAGK superconductivity in MgB2 
by the Akimitsu groupj' triggered enthusiastic studies 
both in experiment and theory to clarifiy the supercon- 
ducting properties and mechanism of the new high-T^ 
superconductor. Given Tc ~ 40K amply exceeding the 
upper bound of Tc ^ 30K which was estimated for con- 
ventional BCS-type superconductors p-* it is a very im- 
portant issue if its superconducting mechanism is con- 
ventional or not. Already Bud'ko et alai has reported 
an isotope effect of a = 0.26 with respect of boron, in- 
dicating the importance of electron-phonon interaction. 
Band calculations have suggested that the diniensionless 
electron-phonon coupling constant A ~ O.TpBS) which 
would be able to give the Tc value of the order of 40K 
if the Coulomb pseudo potential term /z* is very small. 
TunnelinftjCfisiilts indicate the s-wave nature of supercon- 
ductivityQ' B lP All these suggests the BCS-type conven- 
tional superconductivity in MgB2. 

However, a detailed analysis of the specific heat indi- 
cated that the superconducting gap must be anisotropic 
or two-band-like .liH* Non-conventional features were re- 
portecLalso about the power-law temoeijalure dependence 
of i?citif and the penetration depth. tilllj' Theoretically, 
Furukawa pointed out that the Fermi surfaces of two tt- 
bands of MgB2 are close to perfect nesting, which he 
suggests would give rise to antiferromagnetism; the per- 
fect nesting would be realized if holes did not enter lower 
CT-bands, breaking th£.equality of electron and hole num- 
bers in the Tr-bands.liB' In this letter we report that the 
interband polarization function in the MgB2 band even 
as it is greatly enlarged for appropriate wave vectors be- 
tween the two TT-bands and that it can lead to a divergent 
enhancement of the amplitude of Coulomb-origin inter- 
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band electron-pair transfer process. The latter means 
an instability to a two-band-type superconductivity. If 
this mechanism is the dominant one, the type of super- 
conductivity is s-wave-like in each band with each gap 
parameter of an sign opposite between the two 7r-bands; 
the gap in the cr-bands are subsidiary and smaller in 
a clear contrast to the electron-phonon scenario. This 
mechanism allows to qualitatively explain the isotope ef- 
fect when it takes account of the effect of the lattice 
subsidiarily. According to this mechanism Tc can go up 
further through improvement of the nesting. 

Our Hamiltonian is a Hubbard model whose one-particle 
part is the tight-binding bands giving the two 7r-bands 
on the graphene-structure boron layers. Let_us put two 
boron atoms in the unit cell taken by WallaceO at (±0/2-^3, 
0, 0), where a is the distance between the equivalent 
nearest neighbor boron sites in the boron plane. Then 
the tight-binding bands of the 7r-bands are given by 

e ±(k) =±|/i2i(k)|+2t,cos(cfc,) 

= ±i[3 + 2cos(aA:j^) + Acos{Viak.^/2)cos{aky/2)Y/'^ 

+2t,cos{ck,), (1) 

where 

/i2i(k) = t{H-exp(-iai -k) +exp[-i(ai +a2) -k]} (2) 

with ai = (\/3a/2, a/2, 0) and a2 = (0, —a, 0); here t and 
tz are the absolute value of transfer energy between the 
nearest neighbor boron p^-orbitals in the plane and along 
the z-direction, respectively; c is the z-lattice constant. 

The interaction part comes from the on-site Coulomb 
energy Uq at each boron site. The annihilation oper- 
ators at two inequivalent sites are Fourier-transformed 
and here denoted as and bf^^, where a means spin. 
When the one-electron part is diagonalized into two it- 
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bands whose annihilation operators are af^^ and (3}^^, 
the former operators can be hnkcd to the latter through 
the following unitary transformation: 

«fc. = ( exphi(/.(k)](afc^ - /3fc J, 

= (l/^/2) exp[i(/)(k)](afc^ + /5fcJ, (3) 

where (/)(k) is defined by half of the phase of ft.2i(k) in 
Eq. (|^) . If we restrict q appearing in the equation below 
to have only a 2;-component, this restricted part H'j^ of 
the interaction Hamiltonian is rewritten in terms of the 
new operators in a relatively simple form as follows: 

H'^ = {1/N) J2 

ki,k2,q±x,y 















+q^^k2i^ka 


+ (« ^ /?)] 






+qi°ii2i^k,i - 


H (a <-> /?)] 




-q^^k2^ 




- (" (4) 



with 

U = U' = K = L = C/o/2, (5) 

where phase factors cancel each other since 0(k) depends 
only on the in-plane component of k; {a <-> j3) denotes a 
term identical to the preceding one except for exchange 
of a and (3. The term with coupling constant U gives 
rise to the intraband Coulomb scattering, the term with 
U' the interband Coulomb scattering, the term with K 
the interband electron-pair scattering and the term with 
L the interband exchange scattering. This transforma- 
tion is similaii to that performed for the two-chain Hub- 
bard model.ciP The interaction part with q having non- 
z-compon ent is very complicated. 
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Fig. 1. Ladder diagram enhancing the pair transfer process be- 
tween two bands 1 and 2. Straight Unes denote electron Green's 
function. Lateral dashed lines mean interactions with coupling 
constants K and U' . Attached symbols illustrate an example 
among the diagrammatical processes. The lower-case letter de- 
notes a set of wave vector and frequency. 

The interband pair-transfer term ^foXhJK in Eq. (4) is 
known to promote superconductivitytaliJ) If K is larger 
than C/, the simple mean-field calculation gives a finite 
Tc. The amplitude of interband pair-transfer process is 
known to be enhanced due to the higher order processes 
illustrated by the ladder- type diagram in Fig. 1. When 
q in the diagram is chosen parallel to z, the interaction 
denoted by a dotted line is given by simple coupling con- 



stants K or U' and summation of the terms represented 
by this diagram gives the following expression :t3) 

ri2 = A7[i-(i^ + ;7')ni2(q)] (6) 

when the frequency transfer is set to zero for simplicity. 
In this diagram with q || i the phase factor disappears. 

Next, we first evaluate ni2(q) in the case of the perfect 
nesting and then discuss on it in the off-set case. The xy- 
coordinates of the corner points, K (and also H), of the 
hexagonal Brillouin zone are (±27r/-\/3a, ±27r/3a) and 
(0,±47r/3a). At these points the intra-plane part of the 
band ±|/i2i(k)| in Eq. (1) vanishes. So for an appropri- 
ate value of Fermi energy Ei;- in the range between —2tz 
and 2^2, the Fermi surfaces for the ± bands, i.e., electron 
and hole bands, respectively, cross the corner K-H-K line 
at points P± specified by kz = ± cos~^(eF/2tz) = ±fczo, 
repectively. Then, the Fermi surfaces for the electron 
and hole bands in the neighborhood of P+ are given by 
the ± branch of the following equation: 

±|/i2i(k)| = 2tz cos{ckzo) ~ 2tz cos{ckz) 

^ 2ctz sin{ckzo){kz ~kzo). (7) 

As seen in Fig. (2), at least in the neighborhood of point 
P+, where the second equality in Eqs. (7) is justified, 
the electron and hole Fermi surfaces form the mirror im- 
ages with respect of the xy-parallel plane crossing the 
corner line at P+, as seen in Fig. 2 which is shown for 
another purpose. There is another set of approximate 
mirror images around the P_ point. Therefore, if the 
electron Fermi surface just below P_ is translated by 
2kzo in the /c^-direction, it nests with the hole Fermi 
surface just below P_|_. Similarly, the electron Fermi 
surface part just above P+ nests with the hole Fermi 
surface part just below P_ if the former is translated 
by —2kzo- In the neighborhood of the corner line K-H- 
K labeled by (27r/-\/3a,27r/3a) in the fc^^fcy-coordinates, 
expanding /i2i(k) in powers of Px = k^ — 27r/-\/3a and 
Py = ky — 27r/3a, Eq. (7) is approximated as 

± {V3/2)atp\i = 2ctzSm{ckzo)ikz - kzo), (8) 

where p|| = ^Jpx ~^ Py- This shows the Fermi surfaces 
around P+ are approximated by cones as seen in Fig. 2. 
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Fig. 2. Section of model Fermi surfaces with the pajfc^-plane in 
the case of k^o = tt/2c. Cones are considered to correctly repro- 
duce the featurss lef the Fermi surfaces in the neighborhood of 
the cone points,Dt3^ although the peripheral part of the cones 
may not be a good approximation. Hatched part shows the elec- 
tron Fermi surface. The diamond around the center is the hole 
Fermi surface. 

When kzo — tt/2c, one can confirm from Eq. (7) that 
the electron and hole Fermi surfaces take identical forms 
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and perfectly nest with each other by translation of tt/c 
exactly in the fr|amework of the tight-binding bands in 
accord with RefJlj' The interband polarization function 
between the two 7r-bands for zero frequency is defined by 



ni2(Q) 



1 ^ /(£_(k + Q))-/(g+(k)) 

e+(k)-£_(k + Q) ■ 



(9) 



where N is the number of unit cells and f{e±) is the 
Fermi-Dirac distribution function. In order to capture 
gross features, wc approximately evaluate this summa- 
tion in the case of perfect nesting, approximating the 
right-hand side of Eq. (1) by terms linear in p|| and in 
kz — kzo as e±(k) = ±(-\/3/2) aip|| + v±{kz — fczo) with 
v± — 2ctzShi(ckzo) and kzo ~ tt/2c. Then, the Fermi 
surfaces are given as in Fig. 2. The hole Fermi sur- 
face takes the form of two cones sticked at the bottoms 
with total height of tt/c. It centers at kz = 0. The elec- 
tron Fermi surface takes the identical form and centers at 
kz = tt/c in the extended zone scheme. Since at each cor- 
ner integration range is restricted to 120 degrees, we have 
two 360 degree integrals, counting six corners. The up- 
per bound of p|| is set to fc2 = 2(7r/\/3)^/2/a so that the 
total integration area is equal to the area of the reduced 
Brillouin zone. For the wave vector Qo = (0,0, tt/c) of 
perfect nesting, we gets 



ni2(Qo) 



Vcell r, TJ" 



2 vn ^2c' 
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ln(l - 
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}, (10) 



where Vceii = V3a^c/2 is the volume of the unit cell, 
U|| = V3at/2, 7 = 1.78107 and T is the temperature. 
Putting a = 3.0841 and c = 3.5221 and using t - 2.5eV 
and tz ^ 1.5eV,EP 11x2 (Qo) is estimated as 



ni2(Qo) = 0.053 + 0.139 x ln(30152/r)(eV)" 



(11) 



where T is in units of K. Putting T — 40K, we get 
ni2(Qo) = 0.97 (eV)~^. This value is much larger than 
the intra-band polarization function Ili (0) = Vceiiv^/ic^v'^^ 
0.139 (eV)~^ of band i =1 or 2 for wave vector equal 
to zero. For wave vector Q slightly deviating from Qo, 
ni2(Q) decreases in the way that in the argument of 
IninEq. (10) is replaced by -hc|| (Q^ -^Q^) .^^^^(g^ _ 

Qzo)^i where cy and c± are constants. 

In the case of kzo ^ 7r/2c where the sizes of the cones 
become different, only a single cone of two electron Fermi 
cone nests with one of two hole cones for ±Qi = (0, 0, 
±2^20) which is the nesting vector in this case; the other 
combination of cones are not nesting at the same time. 
This situation decreases ni2(Qi) in Eq. (10) mainly in 
the way that the term containing In is divided into two 
terms each with half of the present coefficient; one half is 
not largely changed, keeping the InT singularity, and the 
other diminishes substantially since in the argument 
of In is replaced by _|_ c±{2Tr/c — QziY- 

The magnitude of ni2(Q) evaluated above and its re- 
markable growth around Qz = 2kzo £^re in a fair agree- 
ment with Ease's estimation of II12 (Q) based on FLAP W 
band calculations-EH* However, the prominent peaking 



with decreasing temperature is moderated to a finite 
peak height of 0.265(eV)~"'^. This discrepancy from the 
above obtained divergent behavior is due to the em- 
ployed approximation linearizing £±(k) in p|| and kz — 
kzo- When kzo = ''■/2c, the second order in kz — fczo van- 
ishes in £±(k) and so the linear approximation is not so 
bad and Eq. (10) is valid. However, when kzo 7^ 7''/2c, 
the {kz — fczo)^ term in e±(k) suppresses the InT diver- 
gence, leaving a non-divergent peak around Qz ~ 2fczo. 
In Hase's result there is another peak at Qz — n/c pre- 
sumably due to closeness to perfect nesting. Its height 
is 0.27(eV)~^, slightly higher than the off-set peak. In 
order to decide which is the highest it is necessary to 
carry out numerical work taking account of the depen- 
dences on parameters such as temperature, band filling 
and band parameters. Anyway, ni2(Q) is remarkably 
enhanced around Qi due to partial nesting and around 
Qo due to the closeness to perfect nesting. 

As seen above, Hi2(Q) is very much enhanced around 
Qi and Qo, Eq. (^ is probable to diverge or to be greatly 
enhanced with decreasing temperature, if Uq takes a re- 
alistic value of ^ (2 4)eV since K + U' ~ Uo, from Eq. 
(|^). This means a great enhancement of the amplitude 
of the interband pair transfer process. This leads to a 
superconducting instability due to the two-band mecha- 
nism. The above divergence is the first to appear, among 
other possible divergences, with increasing Uq or with de- 
creasing temperature. 

Tc is determined as the temperature where the pair 
electron scattering diverges in both bands wiiJi. taking 
account of other interaction processes, as in Re^tS Semi- 
quantitatively, however, taking account of the effective 
interband pair scattering coupling constant averaging ri2 
in Eq. (^, by representing it by an effective coupling 
constant K , and also of the intraband Coulomb coupling 
constant U , we obtain the mean- field-like superconduct- 
ing gap equation with two gap parameters Ai and A2 
for the two bands as in Ref.tEP as follows: 



k 



2E,t 



tanh ■ 



E 



k 



2kBT 



1 



A. 



TV ^ ^2E.i 



■ tanh 



E. 



k 



2fcBT' 



(i = 1 or 2) (12) 



where i 



i and S^fc = 



fc 



A?; here we added 



an intraband phonon-mediated attractive coupling con- 
stant Wi in band i as a subsidiary ingredient to consider 
about the isotope effect; is the energy of band i with 
reference to the chemical potential, and is the Debye 
energy. In a simple case where the geometrical means 
of the band half widths of two bands are equal, say to 
D, and the coupling characteristics are identical, the gap 
equation is easily solved to give 



Te = 1.13^13 exp[-l/(A-Ai* 
A = A^i(0)W^i, 

N,{0)iU-K) 



l + N,{Q){U-K)H^) 



(13) 
(14) 

(15) 
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where A^i(O) is the state density per spin of band 1. In 
the present case K is enlarged due to the band nest- 
ing and the ladder process in Fig. 1, /i* becomes small 
and can be even negative and so can be very much 
elevated. It should be noted that the effect of Wi is sub- 
sidiary; evenif = A = 0, = l.lSw^ exp(-l/iVl(0)(ii:- 
^7)) is finite \i K > U . One feature of this model is that 
Al and A2 take the opposite signs although the gap is 
s-wave-like, being constant, in each band. 

The isotope exponent a is given by a = {1 — [^*/(A — 
^*)]2}/2. Applying this and Eq. (13) MgB2 with a = 
0.26 and ujd = TOOK^ we get ^JL* = -0.23 and A = 0.103. 
The latter value looks much smalleiilhan the experimen- 
tally obtained value of A ~ 0.6.Lil'E3' However, the ex- 
perimental A value extracted from the specific heat is an 
upper bound since the mass enhancement due to electron 
correlation is not taken into account. Further, it should 
be noted that our A is for a single 7r-band while the ex- 
perimental A is the sum for all bands. So the above- 
mentioned value of A is not absurd. There is another 
point to note on the isotope effect, i.e., the increase of Tc 
by IK in Mg^*'B2 may have been brought about partly 
through band changes and increasalH) of ni2(Q) due 
expansion of lattice constant a and slight decrease of cl 

In the above consideration only the intraband phonon- 
mediated BCS attraction is taken into account with the 
interband BCS term neglected. The latter is in conflict 
with the K term. This neglect is justified since the elec- 
tron scattering process from one 7r-band to the other by 
the in-plane stretching mode phonons was found to van- 
ish due to symmetry of Bloch wave functions in the close 
neighborhood of the Brillouin zone boundary where im- 
portant carriers exist. 

Any other combination in MgB2, e.g., between a tt- 
and a cr-band, or between two cr-bands, there is no good 
nesting enhancing T\i2 . Further the interband pair trans- 
fer coupling K in these combinations is found to be one 
order of magnitude smaller than Uq. The largest contri- 
bution to K among the cr-bands was found to vanish due 
to symmetry in another tight-binding approximation. So 
in the present framework, the role of the cr-bands are sub- 
sidiary with smaller gap parameters. However, there can 
be an opposite situation, unfortunately to us, where the 
cr-bands are primary due to phonon-mediated attraction 
and the present mechansim of the 7r-bands are not impor- 
tant because the above derived superconducting instabil- 
ity happens to be too weak but it does not look natural 
in view of the jump-wise increase of of MgBc among 
diborides and the large value of ni2 (Q) . Recently Itoada 
also presented a two-band consideration on MgB2 .E3 He 
takes no account of nesting and the two bands in consid- 
eration are different from ours. 

In summary, coupling constant K for the interband 
pair transfer process, coupling constant U' for interband 
Coulomb repulsion etc. in two 7r-bands of MgB2 were 
found to be given by half of the on-site Coulomb en- 
ergy on the boron p^-orbital, specifically in scattering 
processes relevant to the two-band mechanism. In terms 
of them the diagrammatic amplitude for the interband 
pair scattering expressed by ladder diagrams is given by 
Kl[\ - {K + [/')ni2(Q)], where Hi2(Q) is the inter- 



band polarization function between two 7r-bands when 
wave vector Q is parallel to z. Using the tight-binding 
bands for 7r-bands, Hi2(Q) was estimated and found to 
be quite large for Q equal to Qi and Qq which are par- 
allel to z and allow approximate nesting. Therefore, 
-K/[l - [K + C7')ni2(Q)] is very probale to be diver- 
gently enlarged in MgB2. This necessarily brings about 
the two-band-type superconductivity whose Tc can be 
very high. A semi-quantitative gap equation was given 
which takes also account of the BCS interaction in each 
band as a subsidiary ingredient and and the isotope 
exponent were derived. Some characteristic features of 
the gap parameters for various bands were briefly noted. 
This mechanism was claimed to be a very probable candi- 
date for the superconducting mechanism of MgB2. The 
isotope exponent a = 0.26 and the electron-phonon A 
value ~ 0.6, taken by many researchers as supporting 
the electron-phonon scenario, were pointed out not to be 
decisive. The present mechanism gives a guiding princi- 
ple to further enhance T,. that the nesting between two 
TT-bands should be improved. 

The author heartily thanks Dr. I. Hase, Professor J. 
Kondo, Dr. N. Shirakawa and Dr. T. Yanagisawa for 
useful information. 
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